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Abstract 

We study nonlinear stability of spatially homogeneous oscillations in reaction-diffusion sys- 
tems. Assuming absence of unstable linear modes and linear diffusive behavior for the neutral 
phase, we prove that spatially localized perturbations decay algebraically with the diffusive rate 
^-"72 space dimension n. We also compute the leading order term in the asymptotic expansion 
of the solution, and show that it corresponds to a spatially localized modulation of the phase. 
Our approach is based on a normal form transformation in the kinetics ODE which partially 
decouples the phase equation, at the expense of making the whole system quasilinear. Stability 
is then obtained by a global fixed point argument in temporally weighted Sobolev spaces. 
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1 Introduction and main results 



Synchronization of spatially distributed dissipative oscillators has been observed in a wide variety of 
physical systems. We mention synchronization in yeast cell populations [12], fireflies [1], coupled laser 
arrays [24j , and spatially homogeneous oscillations in reaction-diffusion systems such as the Belousov- 
Zhabotinsky reaction [30] and the NO+CO-reaction on a Pt(lOO) surface [29j. Synchronization 
strikes us most when the system size is large or the coupling strength is weak. Both situations relate 
in natural ways to the regime of large Reynolds number in fluid experiments, where one expects 
turbulent, incoherent rather than laminar, synchronized behavior. Still, one finds synchronization 
as a quite common, universal phenomenon, even in very large systems. 

The aim of this article is to elucidate the robustness of spatially homogeneous temporal oscilla- 
tions in spatially extended systems, under most general assumptions, without detailed knowledge 
of internal oscillator dynamics or coupling mechanisms. In fact, quantitative models are very rarely 
available for the systems mentioned above. Instead, we make phenomenological assumptions, related 
to the existence of oscillations and the absence of strongly unstable modes. These assumptions typi- 
cally guarantee asymptotic stability of a spatially homogeneous oscillation in any finite-size system, 
when equipped with compatible (say, Neumann) boundary conditions. The results in this article 
are concerned with infinite-size, reaction-diffusion systems, 

ut = DAu + f{u), u = u{t,x) eM.^ , xGR", t>0, (1.1) 

with positive coupling matrix D G AlArxAf(]R), D = D'^ > 0, and smooth kinetics / G C°°(R^,IR^). 
In this spatially continuous setup, working in the whole space M" is an idealization which corresponds 
to the limit of small coupling matrix and/or large domain size. We will briefly comment on the 
relation between our results in the whole space and the stability of temporal oscillations in finite 
domains, below. 

To be specific, we make the following assumptions on the kinetics / and the coupling matrix D. 

Hypothesis 1.1 (Oscillatory kinetics) We suppose that the ODE ut = f{u) possesses aperiodic 
solution ti*(t) = ti*(t + r) with minimal period T > 0. 

In particular, to avoid trivial situations, we assume that the periodic orbit is not reduced to a 
single equilibrium. As is well-known, this is possible only if > 2, i.e. if the system (II. ip does not 
reduce to a scalar equation. 

In addition to existence we will make a number of assumptions on the Floquet exponents of the 
linearized equation 

Ut = DAu + f'{u4t))u, (1.2) 
which is formally equivalent to the family of ordinary differential equations 

Ut = -k'^Du + f'{u,{t))u, A;GM". (1.3) 

For each fixed k we denote by Fk{t, s) the two-parameter evolution operator associated to the linear 
time-periodic system (II. 3p . so that u{t) = Fk{t,s)u{s) for any t > s. The asymptotic behavior 
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of the solutions of ()1.3p is well characterized by the Floquet multipliers of the system, that is, 
the eigenvalues of the period map Fi.(T,0). We shall rather work with the Floquet exponents 
Ai(A;), . . . , Xwik) € C/icjZ, where m = 2tt/T, which satisfy 



Remark that any Floquet exponent which is simple (that is, of algebraic multiplicity one) depends 
smoothly on the parameter k. Also note that Ai = is always a Floquet exponent for k = 0. We 
refer to the set of Floquet exponents as the Floquet spectrum. 

Hypothesis 1.2 (Marginally stable spectrum) We suppose that the Floquet spectrum in the 
closed half-space {Re A > 0} is minimal. More precisely, we assume that 

(i) The Floquet spectrum in the closed half- space {Re A > 0} is nonempty only for k = 0, in which 
case it consists of a simple Floquet exponent Ai = 0; 

(ii) Near k = 0, the neutral Floquet exponent continues as Xi{k) = —dok'^ + 0(A;'^) for some do > 0. 

We emphasize that these assumptions are satisfied for an open class of reaction-diffusion systems. 
In particular the expansion (ii) with some do S M is a consequence of the simplicity of the Floquet 
exponent Ai = at /c = 0, and of the symmetry k i-^ —k; assuming do > is therefore robust. In 
fact, it is not difficult to show that 



where (•,•) denotes the usual scalar product in M , and f7*(t) is the (unique nontrivial) bounded 
solution of the adjoint equation 



Of course, a necessary condition for Hypothesis [L2] to hold is that u^{t) be a stable periodic solution 
of the ODE ut = f{u), but this assumption alone is not sufficient in general, except if the diffusion 
matrix is a multiple of the identity. Indeed, even if = 2, one can find examples of periodic 
solutions which are asymptotically stable for the ODE dynamics, but become unstable if a suitable 
diffusion is added [211 ESI |23]. One possible scenario, which is usually called phase instability or 
sideband instability, is that the coefficient do be negative, in which case the periodic orbit is unstable 
with respect to long-wavelength perturbations. It may also happen that the Floquet spectrum is 
stable for A; in a neighborhood of the origin, but that there exists an unstable Floquet exponent for 
some k^ ^ 0, and therefore for all /c in a neighborhood of k^. This mechanism is reminiscent of the 
Turing instability for spatially homogeneous equilibria. In Section [7| below, we give an example of 
a simple 2-species system which exhibits both kinds of instabilities depending on the choice of the 
parameters. 

In order to state our results, we introduce a function space which measures both the spatial 
localization and the amplitude of the perturbation to our spatially homogeneous profile . We will 
consider initial perturbations in the space of functions X = L^{W-) n L°°(R"), with target space 
M^, equipped with the norm 




J^{U,{t),DuUt))dt 
/o^([7*(t),<(t))dt 



(1.4) 



Ut = f'iu,{t)fU. 



(1.5) 
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where "sup" here refers to the essential supremum. We will measure decay in the space L°°(M"'). 
Our first result is: 

Theorem 1 Consider a reaction- diffusion system (jl.ip on with oscillatory kinetics, Hypothe- 
sis and marginally stable spectrum, Hypothesis Then there are C,5 > such that for any 
initial data u(0,x) = u*(to) + Vo{x) with to €^ M arbitrary and \\vq\\x < S, there exists a unique, 
smooth global solution u{t,x) of (jl.f p for t >0. Moreover u[t,x) converges to the periodic solution 
in the sense that 



sup 



u{t,x)-uSo + t) <-^^^^, forallt>0. (1.6) 



We emphasize that the perturbations we consider are localized in space, and therefore do not 
alter the overall phase to of the periodic solution. We refer however to Section [7| for a discussion of 
possible stability results in more general situations. It is not difficult to verify that the decay rate 
in (jl.6p is optimal. In fact, under the assumptions of Theorem [H one can even compute the leading 
term in the asymptotic expansion of the perturbation as t ^ +cxd. Let 



where do > is defined in (II. 4p . Our second result is: 

Theorem 2 Under the assumptions of TheoremUl the solution u{t,x) of (II. ip can be decomposed 
as 

u{t,x) = u^{to + t) + u'^{to-^t)a{t,x) f3{t,x) , x G M" , t>0, (1.8) 
where a : R+ x ^ M and /3 : IR+ x ^ satisfy 

||/3(t,-)||Li+(l + t)"/'ll/3(t,-)||L- - 0, (1.9) 

I— ► + 00 

||t"/2a(t,a;\/t) -a=.G||iinLoo > 0, (1.10) 



for some a* G M. In addition. 



in (U^{to),K{to)) 



where U^{t) is the bounded solution of the adjoint equation (jl.Sp . 
In other words, the solution u(t, x) of (jl.ip satisfies 



u{t, x) = u.{to + t)+ Kito + t) e-l-l /(^'^o*) + o{t--l^) 



To leading order, the effect of the perturbation is thus a spatially localized modulation of the phase 
of the periodic solution. As is clear from the proof, the left-hand side of (jl.9p decays at least like 
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as t ^ oo, for some 7 > 0. However, to specify a convergence rate in (jl.lOp . it is necessary 
to restrict ourselves to more localized perturbations. For instance, if we assume in addition that 
(1 + |2;|)?Jo G L^{W^), then we can prove that the left-hand side of (jl.lOp is 0{t~^^'^) as t — > 00. 

Under Hypothesis [L2l if we consider system (jl.ip in a large bounded domain (say x G Q/e where 
0, C M" is bounded) with Neumann boundary conditions, the perturbations of the periodic solution 
u^{t) decay exponentially [13]: if ||vo||l°= ^ ^ for some small 6 > 0, then ||ti(t)||/^oo < C||?;o||l°° e~^* 
for some /x > 0. However, the relaxation rate and the size of admissible perturbations 6 both 
depend on e, with typical scalings fi,S = O(e^) predicted by the spectral gap of the Laplacian on 
the domain Q/e. This gap vanishes in the limit e ^ 0, and Theorem [1] shows that exponential decay 
is replaced by diffusive decay. Nevertheless, we expect our results to give an accurate description of 
the intermediate asymptotics for large bounded domains, if the initial perturbations are sufficiently 
localized. 

The type of diffusive decay that we establish in Theorems [1] and [2] has been observed in many 
other contexts. For instance, localized perturbations of spatially periodic, stationary patterns in the 
Ginzburg-Landau or Swift-Hohenberg equation exhibit a similar diffusive behavior [2ll5|[6l [T5l[26| [28]. 
At a technical level, the approach in [21 [261 HH] is based on renormalization group theory, see for 
instance |3|. Roughly speaking, the method relies on the fact that the time-T map for the evolution 
of the perturbations becomes a contraction in a space of localized functions when composed with 
an appropriate rescaling, except for a neutral direction which specifies the profile of the self-similar 
solution describing the leading order asymptotics. A nice feature of renormalization theory is that 
it allows to determine easily which terms in the nonlinearity are "relevant" (that is, potentially 
dangerous) for the stability analysis. For example, if we consider the nonlinear heat equation 
Ut = Au + \u\P in M", with small and localized initial data, it is well-known that the nonlinearity 
I up will not infiuence the decay predicted by the linear evolution Hp > 1 + 2/n, whereas instabilities 
and even blow-up phenomena can occur if p < 1 + 2/n [10^ [6], In particular, quadratic terms (which 
arise naturally in the Taylor expansion of any smooth function) are "irrelevant" if n > 3 and 
"relevant" if n = 1,2. For this reason, diffusive stability is often easier to establish in high space 
dimensions, when diffusion is strong enough to control all possible nonlinear terms, whereas serious 
problems can occur in low dimensions. This is the case in particular in the stability analysis of one- 
dimensional spatially periodic patterns [261 19]. where a key step of the proof is to show that relevant 
"self-coupling" terms actually do not occur in the evolution equation for the neutral translational 
mode. 

As one may expect from the discussion above. Theorem [1] is rather easy to prove when n > 3. For 
completeness, we first settle this case in Section [2] and then focus on the more interesting situation 
where n = 1 or 2. Here the idea is to construct a normal form transformation for the ODE dynamics 
which removes all "relevant" terms in the nonlinear PDF satisfied by the perturbation. In Section [S] 
we show that this is possible, at the expense of transforming the semilinear equation (|l.ip into a 
quasilinear parabolic system. The next important step is to obtain optimal decay estimates for the 
solutions of the linearized perturbation equation, including maximal regularity estimates, using the 
spectral assumptions in Hypothesis II. 2[ Since the perturbation equation is translation invariant in 
space and periodic in time, such bounds are relatively straightforward to obtain via Fourier analysis, 
see Sectional Using these linear estimates, we give in Section[5]a proof of Theorem[T]which is valid for 
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n < 3, hence covering the missing cases n = 1,2. Instead of renormalization group theory, we prefer 
using a global fixed point argument in temporally weighted spaces, as in [5l[6l[T5]. After stability 
has been established, a rather classical procedure, which is recalled in Section [6l allows to derive 
the first-order asymptotics and to prove Theorem [2] at least for n < 3 (the higher dimensional case 
is again easier, and left to the reader). In the final Section [71 we illustrate our spectral assumptions 
on a simple, explicit example, and we conclude with a short discussion including possible extensions 
of our results. 

Acknowledgments. The authors thank Sylvie Monniaux for useful discussions on maximal regu- 
larity in parabolic equations. A.S. would like to thank the Universite de Franche-Comte for generous 
support and hospitality during his stay, where part of this project was carried out. A.S also ac- 
knowledges partial support by the NSF through grant DMS-0504271. 

2 Stability in high dimensions 

In this section we explore a straightforward and somewhat naive approach to the stability of the 
periodic orbit u,,(t) as a solution of the reaction-diffusion system (jl.ip . This method gives a simple 
proof of Theorem [1] in the high-dimensional case n > 3, the main ingredient of which is an L^-L'^ 
estimate for the linearized evolution operator which will be established in Section HI Without loss 
of generality, we assume from now on that the parameter to in Theorems [T] and [2] is equal to zero 
(this is just an appropriate choice of the origin of time). 

Consider a solution u{t,x) = n*(t) + v{t,x) of (|l.ip . The perturbation v satisfies the equation 

vt = DAv + f'{u^{t))v + N{u4t),v) , (2.1) 

where 

N{u,{t),v) = f{u,{t)+v)-fiu4t))-f'{u4t))v. 

The Cauchy problem for the semi-linear parabolic system (j2.ip is locally well-posed in the space 
X = Li(M") n L°°{W), see e.g. [I3l[18]. More precisely, for any vq E X, there exists a time f > 
(depending only on ||fo||j>c) such that (|2.ip has a unique (mild) solution v S C^{[0,T], L^{W^)) n 
CftO((0,r],L°°(M")) satisfying v{0) = vq. 

Remark 2.1 Due to parabolic regularization, the solution v{t,x) of (j2.ip is smooth for t > 0. For 
instance, there exists C > such that \\v{t)\\u2 < Ct~^\\vo\\x for all t G (0, T]. Therefore, in the 
proof of Theorem {1\ we can assume without loss of generality that the initial perturbation is small 
in the space X n H"^ (W) . 

To investigate the long-time behavior of the solutions of (j2.ip . we consider the corresponding 
integral equation 

v{t) = T{t,0)vo + J^{t,s)N{u^{s),v{s))ds, (2.2) 
Jo 

where J^{t,s) is the two-parameter semigroup associated to the linearized equation (jl.2p . Due to 
our spectral assumptions (Hypothesis II. 2p . the operator J^{t,s) satisfies the same L^-L'i estimates 
as the heat semigroup e^*~*^^. More precisely, anticipating the results of Section [H we have: 
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Proposition 2.2 {L^—L^ estimates) There exists a positive constant C such that, for all t > s 
and all 1 < p < q < oo, we have 

(J 

\\^{t,s)v\\L^i^M.r^) < nn.i) \Mlp{R")- (2-3) 

(t-s)2*'P q> 

Proof. The proof follows exactly the same lines as in Propositions 14. 1^ 14.31 and 14. 4[ ■ 

By construction, the nonlinearity N{u^:,v) in (|2.ip is at least quadratic in f in a neighborhood 
of the origin. More precisely, there exists a nondecreasing function K : R_|- such that, for all 

t£[0,T], 

\N{u^{t),v)\ < K{R)\v\'^ whenever <R. 

As was mentioned in the introduction, if the space dimension n is greater or equal to 3, the diffusive 
effect described in (12. 3p is strong enough to kill the potential instabilities due to the nonlinearity. In 
that case, nonlinear stability can therefore be established by a classical argument, which we briefly 
reproduce here for the reader's convenience. 

Proof of Theorem [U (n > 3). Fix vq e X = ^^(R") n L°°(M"), and let v G C^{[0,T^), L^{W)) D 
C°((0,r*),L°°(M'^)) be the maximal solution of I^Jij with initial data vq. For t G [0,r*) we denote 

= sup \\v{s)\\li+ sup + . 

0<s<t 0<s<t 

Using the integral equation (j2.2p and the linear estimates (j2.3p . we easily find 

\\vit)\\Li < mt,0)vo\\Li + / mt,s)N{u,{s),v{s))\\Lids 

Jo 

<C\\vo\\L^+CK{^{t)) I \\v{s)\\li\\v{s)\\l^ ds 

Jo 

<C\\vo\\L^+CK{cPitmtf __L_ds. 

Similarly, if < t < 1, we have 

(1 + (OIIloo < Cllt-ollioo + C7 r ||iVK(s), v{s))\\l^ ds 

Jo 

<C\\vo\\L^+CK{4>{t))4>{tf [ ^ ds, 

Jo + 

while for t > 1 we can bound 

(1 + tr/^vit)h^ < c\\voh^ + c{i + tr'^ j'^' ^t-sT^ 

+ C{l + tr'^ f \\N{u,{s),v{s))\\Lo^ds 

Jt/2 
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Now, since n > 3, we have /o°°(l + s) ""/^ ds < oo and we see that there exist positive constants 
Ci,C2 (independent of T^) such that 

Ht) < <^i||^olU + C2K{(l){t))(P{tf, for ah t £ [0,T,) . (2.4) 

So if we further assume that the initial perturbation vq £ X is smah enough so that 

2Ci\\vo\\x < 1, and 4CiC2K{l)\\vo\\x < 1, 

then it follows from (j2.4p that (/)(t) < 2Ci||wo||x < 1 for all t G [0, T*). Since [0,T^) is the maximal 
existence interval, this bound implies that T^, = +oo and that the solution of (j2.ip satisfies 

sup||t;(t)|Ui +sup(l + t)"/2||^(t)||^^ < 2C7i||t;o||x. 

t>0 t>0 

This concludes the proof of Theorem [1] in the high-dimensional case n > 3. □ 

3 Reduction to a normal form 

In low space dimensions the argument presented in the previous section fails, and we must therefore 
have a closer look at the structure of the perturbation equation. The idea is to introduce a normal 
form transformation which simplifies the ODE dynamics in a neighborhood of the periodic orbit 
ti*. Applying this transformation to the reaction-diffusion equation (jl.ip . we obtain a quasilinear 
parabolic system which will be the starting point of our stability analysis in Sections U] and [5l 

We thus consider the ordinary differential equation 

ut = f{u), ueR^, (3.1) 

with smooth nonlinearity / G C°°(M^,M^), and we assume the existence of a time-periodic solution 
w*(t) = u^{t + T) with minimal period T = Itt/lo > 0. As in Hypothesis [L2l we suppose that is 
linearly asymptotically stable, in the sense that the Floquet exponents Ai, . . . , Xn are all contained 
in the open left half-plane, except for Ai = (which is therefore algebraically simple). We shall show 
that the dynamics of ()3.ip near the periodic orbit ti* is conjugate to the dynamics of the following 
normal form 

9t = Lu, vt = g{B,v), 0G5^^M/27rZ, veB^cR^'^, (3.2) 

where denotes the open ball of radius e > centered at the origin in M^^^. Here the vector field 
g has the expansion 

g{9,i) = L{e)i + g2{9,v)[v,v], (3.3) 

where L{9) is a real (A^ — 1) x (A — 1) matrix depending smoothly on 6, and (72(^5 v) is a symmetric 
bilinear form on M^~^ depending smoothly on 6, v. In particular g{9, 0) = 0, hence (|3.2p has a 
trivial solution 6(t) = ujt, v{t) = which will correspond to the periodic solution n^,(t) of (j3.ip . By 
construction, the Floquet exponents A2, • • • , Aat of the time-periodic linear operator L{ujt) are all 
contained in the open left half-plane. 

In what follows we denote by $(t) the flow of (|3.ip in a neighborhood of n^,, and by $nf(i) the 
flow of (j3.2p in a neighborhood of x {0}. These local flows are defined at least for t > 0. 
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Lemma 3.1 (Normal form) Assume that the periodic solution n* is linearly asymptotically stable. 
Then there exist e > e' > and a smooth diffeomorphism ^ from the solid torus x to a tubular 
neighborhood of the periodic orbit such that the local flow in x M^"-"^ defined on x B^i by 

^,,^{t) = ^'"^ o o , t > , 

is the flow induced by an ODE of the form (j3.2p , (j3.3p . 

Proof. Since the periodic orbit is linearly asymptotically stable, we can find a tubular neigh- 
borhood which is smoothly foliated by strong stable fibers. Straightening out these fibers gives 
the desired representation of the fiow. For completeness we construct this straightening change of 
coordinates in detail. 

We start with the linearized equation at the periodic orbit, ut = f'{u^:{t))u, which possesses a 
linear invariant smooth foliation: if we parametrize the periodic orbit u^:{t) using 9 = ujt ^ , 
Floquet theory gives smooth families of complementary subspaces E^^{9) and E^{9), such that 
dim.E^^[9) = N — 1 and E^{0) = span (u^^/o;)). The linearized evolution leaves these subspaces 
invariant: u{t) £ E^'^^O) implies u(i + r) G E^'^^O + lot), and the same holds for E^{9). In particular, 
the family {E^^(9)}Qfzs^ forms a smooth normal bundle to the periodic orbit (which is an orientable 
manifold in ambient Euclidean space), and such a bundle is necessarily trivial. Thus, we can find 
smooth coordinates {9, v) €z x R^"-*^ and a smooth map ^-0 : 5"^ X M^-i ^ such that 
^o(^,0) = u,{9/uj) and '^{9,R^-^) = u^{9 / u) + E^\9) for all 9 G . 

On the other hand, for each 9 £ S^, the strong stable manifold W'^^{9) of the nonlinear system 
(1311) is the graph of a local map he : ^'"(6') E''{9), with he{0) = and h'g{0) = 0. The strong 
stable manifolds depend smoothly on the base point. In other words, hg depends smoothly on 9, so 
that the map 

u,i9/Lu) + v'' ^ ^i{u,{9/iu)+v'') := u,{9/u;) + v'' + he{v'') , 

defines a smooth difi'eomorphism in a tubular neighborhood U of the periodic solution. Thus, if 
e > is sufficiently small, the map ^ := o : x ^ U is also a smooth diffeomorphism 
onto its image, and ^{9,B^) C ^""{9) for all 9 G S^. Since 

<^{t){W^%9))nU C W''''{9 + ujt), (3.4) 

we deduce that, i( 9 £ and v G B^i for some small e', then (<&(t) o '^)(9,v) belongs to the image 
of ^ for ah t > 0, and 

<^^i{t){9,d) := {^f-^ o<^{t)oqf){9,v) = {9 + ut,v), 

for some v £ B^. This immediately implies the trivial form 9t = to for the evolution equation 
associated to $nf- Moreover, by construction, <I>nf(0,O) = {9 + u)t,0) for all t > 0, hence the 
transverse variable v evolves according to an ODE of the form (j3.2p . where the vector field satisfies 
g{9, 0) = and can therefore be expanded as in (|3.3p . ■ 

We conclude this section with the transformation of the full reaction-diffusion system. The 
pointwise change of coordinates u = ^(v) = ^{9,v) yields 
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which can be expanded into 

vt = 'i>'{v)-^D'i>'{v)Av + '^'{vy^D'^"{v)[Vv,Vv]+ fniiv). (3.5) 

We are interested in the stability of the periodic orbit v^{t) = (wt, 0)^, and therefore we set v = 
v*{t) + w{t,x), so that w solves 

wt = '^Xv,+w)-^D^\v^ + w)Aw + ^\v,+w)-^D^'\v,+w)[Vw,Vw]+f^i{v, + w), (3.6) 

where now f^f{v) = {0,g{v))'^. In what follows, the first component of the vector w will play a 
distinguished role, as is clear from the expression of Thus we shall often write w = {wo,W}^)'^, 
with Wo G M and G K^"^ 



4 Linear evolution estimates 



We consider the linearization of (j3.6p at w = 0, which reads 

Wt = ^'{v,)-^D^'{v,)Aw + f;,i{v,)w . (4.1) 
To simplify the notations, we define 



A{t) = "^'(v^t)) , and B{t) 




L(wt), 



where L{9) is the (A^— 1) x (A^— 1) matrix which appears in (|3.3p . Note that A{t), B{t) are T-periodic 
N X N matrices, and that A{t) is invertible for all t. The linearization ()4.ip then becomes 

Wt = A{t)-'^DA{t)Aw + B{t)w, t£M. (4.2) 

By Fourier duality this system is equivalent to the family of ODEs 

Wt = -k'^A{t)-^DA{t)w + B{t)w , tGM, G . (4.3) 

Since (j4.3p is related to (jl.3p by the T-periodic linear transformation u = 'if'{v^(t))w, it is clear that 
the Floquet spectrum of (j4.3p is identical to that of (jl.3p and therefore satisfies Hypothesis II. 2[ 
Let M{t,s]k) denote the evolution operator defined by (|4.3p . so that any solution of (|4.3p satisfies 
w{t) = M{t, s; k)w{s) for t > s. As w = {wo,Wh)'^ G M x it is natural to decompose the 

matrix M in blocks as follows: 



M{t,s;k) 



fMoo{t,s;k) Moh{t,s;k)\ 
\Mi,o{t,s;k) Mhh{t,s;k)J ' 



where Mqo, Moh, Mho, Af^h are matrices of size 1 x 1, 1 x (iV - 1), (A^ - 1) x 1, (iV - 1) x (A^ - 1), 
respectively. The main result of this section is the following pointwise estimate on M{t, s; k): 
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Proposition 4.1 (Pointwise estimates) There exist constants C,d > such that, for all t > s 
and all k G M", one has 

|Moo(t,s;A;)| <Ce-'^^'(*-^), (4.4) 

|Moh(t, s;k)\ + |Mho(t, s;k)\< f ^-'^^\t-s) ^ ^^^^^ 

1 + t — s 

l^hh(t,s;A:)| < e-'^^'(*--) , (4.6) 

(1 + t - s)^ 

where the norms on the left-hand side are arbitrary, k-independent matrix norms. 

Proof. Since the coefficients in (j4.3p are T-periodic, we have M{t + T, s + T; k) = M{t, s; k) for 
ah s G M and ah k € M". As a consequence, if t > s and if ti,T2 G [0, T) are such that t — ti and 
s + r2 are integer multiples of T, we have the identity 

M{t, s; k) = M{t, t-n;k) M{k)"' M{s + T2, s;k), A; G R" , (4.7) 

where M{k) = M(T, 0; k) and m G N is such that t — s = ti + mT + T2. To prove Proposition 14.11 
it is therefore sufficient to estimate M{k)"^ and M(t, s; k) for Q <t — s <T. 

Step 1: Estimates on M{k)'^. 

The general strategy is to distinguish between various parameters regimes. For small and interme- 
diate values of k, we essentially exploit Hypothesis 11.21 while for large k it is sufficient to use the 
parabolicity of (14. 3p . 

Small k: We solve the ODE (j4.3p perturbatively for t G [0,T] and obtain 

w{t) = U{t,0)w{0) -k^ [ U{t,s)A{s)-^DA{s)U{s,0)w{0)ds + O{k^), 



JO 

where U{t,s) = M{t,s;0) is the evolution operator associated to the equation wt = B{t)w. In 
particular, setting t = T, we find 

M{k) = U{T,0) f id f {A{t)U{t,0))-^D{A{t)U{t,0))dt + O{k'^) 







1 0\ 1 - doTk"^ + 0{k^) 0{k^) 
v] [ 0{k^) id +0{k^) 



(4.8) 



where V is the (A^ — 1) x — 1) Floquet matrix associated to the T-periodic linear operator L{ujt), 
and 



do = ^ r eJ{A{t)U{t,0))-^D{A{t)U{t,0))eidt = ^ C elA{t)-^DA{t)eidt. 
^ Jo J- JO 



(4.9) 



Here ei = (1,0)^ is the first vector of the canonical basis in M^. As was already observed, all 
eigenvalues of V are contained in the disk {2: G C | |z| < e"'^} for some > 0, and it follows from (j4.8p 
that M{k) has exactly A^ — 1 eigenvalues in this disk if k is sufficiently small. The remaining Floquet 
multiplier has the expansion 1 — doTk'^ + 0(A;^), in agreement with Hvpothesis 11.21 Incidentally, 
we observe that (j4.9p is identical to (|1.4p . Indeed, since M*(t) = '^{v^{t)) = ^(ujtei), we have 
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u'^(t) = LJ^'{v^:{t))ei = LL!A{t)ei, and it is also straightforward to verify that the bounded solution 
of the adjoint equation (II. 5p is U^,{t) = {A{t)~^)'^ei. Thus ()4.9p can be written as 

do = —p^ I UAt) DuAt)at = , 

and Hypothesis 11.21 guarantees that do > 0. 

For k £ M" sufficiently small, let P{k) denote the spectral projection onto the one-dimensional 
eigenspace of M{k) corresponding to the neutral Floquet exponent Xi{k) = —dok"^ + 0(/c^). From 
it is easy to verify that P{k) has the following form 



P{k) 



1 I k^b 



2iT 



1 + k^T^ I ^2^ ^4^5^ 



where a{k), b{k) are (A^ — l)-dimensional vectors with a{k), b{k) = 0(1) as A; ^ 0. By construction, 
we have for any m G N*: 

Mik)"" = M{k)"'P{k) + M(A:)'"(id - P{k)) = e'"'^^i(*^)p(A;) + (M(A:)(id - P(A;))r • (4.10) 

Since Ai(A;) < —dik"^ for small A; if < di < do, and since the spectral radius of M(/c)(id — P{k)) 
is smaller than e~'^ , we conclude that 

\Mikr\ < Ce-'^^^'-^ (^l, +0{e—) < C e''^^''-^ (^^^^^_, ^^^^''^ , 

for all m G N* if < kq <ti 1. Here the matrix norm | • | is applied separately to each of the four 
blocks of M{k)"^, and for convenience the four upper bounds are collected in a 2 x 2 matrix. 

Large k: In this parameter regime, it is more convenient to set v{t) = A[t)w{t) and to solve the 
v-equation corresponding to (j4.3p . namely 



vt = -k^Dv + C{t)v, where C{t) = ^ {t)A{t)-^ + A{t)B{t)A{t)~^ . (4.11) 
The matrix C{t) is T-periodic, hence uniformly bounded. A standard energy estimate yields 

^^ll^Wlli = -k\v{t),Dv{t)) + {v{t),C{t)v{t)) < -d2k^\\v{t)\\l + K\\v{t)\\l, 

where (i2 > is the smallest eigenvalue of the (symmetric and positive) matrix D, and K = 
sup^gjQji] ||C(t)||2. Thus any solution of (|4.1ip satisfies ||v(i)||2 < e(~'^2fc^+ii')i||^(^o)||2, and returning 
to the ty-equation we obtain 

lk(t)||2 < Ce(-'^2^''+^)*||«;(0)||2, t>0, (4.12) 

for some C > 0. In particular, if we choose t = mT and if we assume that > ki with ki = 
i^Kjd'if'l'^ , we arrive at 

\M{kT\ < Ce-'^^'^'"^, where = ^ ■ 
12 



Intermediate k: By Hypothesis 11.21 if kq < \k\ < ki, the spectrum of M{k) is entirely contained 
in the disk {z € C | |2:| < e~'^} for some /x > 0. Moreover, the resolvent matrix {z — M{k))~^ is 
uniformly bounded for all z on the circle {\z\ = e~^} and ah k in the annulus kq ^ 1^1 ^ f^i- If 
< d4 < fx/inlT), it fohows that 

\M{k)"'\ < Ce-"'*'^"'"^, mEN, 

where the constant C is independent of k. 

Summarizing the results obtained so far, we proved that there exist C > and d > such that 



\M{kr\ < Ce-'^^ ™^ ( ' _i ' (4-13) 



for all k £ M" and all m G N*. This is the particular case of (|4.4p - (j4.6p when s = and t = mT. 

Step 2: Estimates on M{t, s; k) for < t - s < 2T. 
Our goal is to show that 

\M{t,s;k)\ < Ce-'^'=^(*-) (^^,^^'_ k\t-s)\^ ^^^^^^ 

for some C > and ti > 0. Note that (iniD imphes (gSD-dlS]) if < t - s < 2r. 

In the case where k'^{t — s) is small, say /c^(t — s) < K2 ^ 1; we can solve ()4.3p perturbatively 
and obtain as in (j4.8p 

M(t,s;fc) = (id+0(A;2(t-s))), (4.15) 

from which (j4.14p follows (for any fixed d > 0). li k"^ {t — s) > K2, then using (j4.12p we immediately 
find 

\M{t,s-k)\ < C7e(-^2'='+^)(*-^) < Ce^^^e-'^^fc^Ci-s)^ 
which implies (j4.14p with d = d2- 

It is now straightforward to conclude the proof of Proposition 14. li In view of (j4.14p . it remains 
to prove ()4.4p - ()4.6p for t — s > 2T. We decompose t — s = ti + mT + T2 with ri, r2 G [0, T) and 
m G N*, and we factorize M(t, s; /c) as in (fiTT]) . Using (|4.13p . ()4.14p . we find 



which is the desired estimate. 

Remark 4.2 For a// A; G W" and all s G M, we /ia?;e 



1 /tVA / 1 (mr)-i\ / 1 k'^T2 



1 (l + t-s)-i\ 

1 + t-s)-^ (l + t-s)-^ ' 



hm M(t,s;fct-i/2) = e-'^°'=M ^ ° 

t^+oo ^ ' ' ^ loo 
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where do is as in Hypothesis This follows from the proof of Proposition \4-l\ and in particular 
from (gZD, Km . 

Proposition 4.3 (Estimates for the derivatives) For any a G N" there exists C > such that, 
for all t > s and all k G M", 

\d^M{t,s;k)\ < C7(t-s)l°l/2e-'^'='(*-«). (4.16) 
The same estimate holds for Mqq, (1 + t — s)Moh, (1 + t — s)Mho, and (1 + t — s)2Mhh- 

Proof. It is clear from (|4.3p that M(t, s; k) depends on the parameter k G only through the 
scalar quantity p = k^. In this proof, we set M(t, s; k) = M{t, s; k'^) and we consider the derivatives 
of M{t, s;p) with respect to p. Our goal is to prove the estimate 

\dlM{t,s;p)\ < C,{t-sye-''P^'-'\ jGN, (4.17) 

which implies immediately (j4.16p . Differentiating (j4.3p with respect to p = k'^, we find 

{dpw)t = -V{t)w + {B{t)-pV{t)){dpw), 

where V(t) = A(t)-'^ D A(t) . Since dpM{s,s;p) = 0, we deduce that 

dpM{t,s;p) = - M{t,ti;p)V{ti)M{ti,s;p) dti . 
Iterating this procedure, we obtain for any j G N the representation formula 



t rt 



dlM{t,s;p) = {-iyj\ / .../ M{t,ti;p)V{ti)x 

X M{ti,t2;p)V{t2) . . . M{tj^i,tj;p)V{tj)M{tj,s;p) dtj . . . dti . 

Using now the pointwise estimates established in Proposition 14.11 we easily obtain (|4.17p . Similar 
bounds can be proved for Mqq, (1 + t — s)Moh; {1 + t — s)MhO) and {1 + t — s)2Mhh (we omit the 
details). ■ 

We now convert our estimates on the Fourier multipliers M(i, s; k) into bounds on the linear 
evolution equation (j4.2p in various spaces. The two-parameter semigroup M{t,s) associated to 
()4.2p is defined using Fourier transform by the relation 

{Mit^)v){k) = M{t,s;k)v{k), keW. (4.18) 

The following proposition contains the main estimates on M.{t, s) which will be used in the nonlinear 
stability proof. 

Proposition 4.4 (L^-L'^ estimates) There exist a positive constant C such that, for allt>s and 
all 1 < p < q < oo, one has 

(J 

\\M{t,s)v\\Li(9.^) < \\v\\lp{r^)- (4.19) 

(t_s)2l.p q> 



The same estimate holds for A^oO; (1 + t — s)A^oh) (1 + t — s)A1hO; '^'^^ {1 + t — s)"^]^ 



hh- 
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Proof. By construction M{t,s) is the convolution operator with the function x M{t, s; x), 
which is just the inverse Fourier transform of k i— > M(t,s;k). Thus using the pointwise estimate 
(j4.4p . we easily obtain 

||A^(t,S;-)||L-{M") < C||M(t,g;-)||Li(R") < (^_^^n/2 " 

To estimate the norm of A4{t,s;-), we use Sobolev embeddings. Let m G N be the smallest 
integer such that m > n/2. Using the estimates of Proposition 14.31 together with Holder's inequality 
and Parseval's identity, we find 

f / \x\'^ \ —m/2 / IxP \m/2 

\\M{t,s;-)\\mur^) = 1 + ^ 1 + ^ \M{t,s;x)\dx 



<C{t-sr/'ll^^ ^{t-sr^-^\d^M{t,s;k)fdk\ <C. 



Summarizing, we have shown that 



C 

and (j4.19p follows by Young's inequality. The estimates for TWqO) (1 + t — s)7Woh5 (1 + i — s)A4ho, 
and (1 + f — s)^A^hh are proved in exactly the same way. ■ 

Remark 4.5 Under the assumptions of Proposition we also have 

C 

||(9^A1(t,s)?;||i,(Rn) < „ ^ lbllLP(E") , aGN", 

and the same estimates hold for , + t — s)A^oh, (1 + ^ — s)-A^hO; and (1 + i — s)^A^hh- This is 
obvious in view of Proposition \4T1\ since the operator d"Ai{t,s) is the Fourier multiplier associated 
to the function {ik)°'M{t, s; k). 

Remark 4.6 For later use, we also observe that, if t > and < s <t, then 

C s 

\\iMoQit,s) - Mooit,0))v\\Li(m") < T Iblk^K") ' (4-20) 

{t - i' ^ 

for 1 < p < q < oo. If t/2 < s < t, this bound follows immediately from ()4.19p and the triangle 
inequality. If < s < t/2, we observe that 

Moo{s,0;k) = 1 + k'^sRo{s,k) , Mho{s,0; k) = k'^sRh{s,k), 

where Ro{s, k), R\i{s, k) are uniformly bounded for s > and k £ R", see (j4.15p . Since 

Moo(i, s; k) - Moo{t, 0; k) = Moo{t, s; k){l - Moo(s, 0; k)) - Moh(t, s; k)Mi,o{s, 0; k) , 
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we obtain the pointwise bound 



\Moo{t,s;k)-Moo{t,0;k)\ < Ck'^s e-'"^^'-'^ < C^q-^ 

t — s 

which allows to establish the U'-L'^ estimate (j4.20p using the same arguments as in the proof of 
Proposition \4-4\ 



Finally, to control the quasilinear terms in the perturbation equation we will use maximal regu- 
larity properties of the evolution semigroup M{t, s). 

Proposition 4.7 (Maximal regularity of type U) For any r G (l,+oo) and any T > 0, there 
exists C > such that the following holds. If v £ L''((0, T), -L^(M")) and if w satisfies 

ft 



w{t) = [ M{t, s)v{s) ds , te[0,f], (4.21) 
Jo 



then 

[ \\Aw{t)\\l2 dt < C [ \\v{t)\\l2dt. (4.22) 
Jo Jo 

Proof. For any i € M the generator £(t) = A{t)-^DA{t)A + B{t) is an elliptic operator in L'^(W^) 
with (time-independent) domain i7^(M"). Moreover C{t), considered as a bounded operator from 
H'^(W^) into L^(M"), is a smooth function of t. Thus estimate (j4.22p is a particular case of the 
results established in |T4l[l9]. ■ 



We also state a corollary which will be useful in the next section. 
Corollary 4.8 Fix r £ (l,oo) and q G M. There exists C > such that, if 

wit) 

then 



[ M{t, s)v{s) ds , t>T, 

/ {i+tr\\Awit)\\i,dt < c {i+tr\\vit)\\i2dt. 

Jt Jo 



Proof. Fix G N* := N \ {0}. If t G [kT, {k + l)r], where T > is the period of C{t), we can 
write ^ ^ ^ 

w{t) = / M{t,s)v{s)ds - / M{t,s)v{s)ds . 

J{k-l)T J(k-1)T 

Since M{t + T, s + T) = M{t, s), and since 0<t-{k-l)T < 2T, we can use Proposition O (with 
T = 2T) to control both terms in the right-hand side. We obtain 

\\Awit)\\l,dt < C / \\v{t)\\l,dt, 

kT J(k-l)T 

for some C > independent of k. Multiplying both sides by (kT)" ~ (1 + t)" and summing over 
A; G N*, we obtain the desired result. ■ 
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5 Nonlinear stability in low dimensions 



This section is devoted to the proof of Theorem [T] in the case n < 3. Smoothness of the change 
of coordinates ^ imphes that it is sufficient to prove the decay estimate (jl.6p for the transformed 
equation (j3.6p . with smahness assumptions on the initial data wq. To simphfy the notations, we 
rewrite (|3.6p in the compact form 



wt = C{t)w + F{t, w, Aw) + G{t, w, Vw) + H{t, w) , (5.1) 

where C{t) = A(t)~^DA{t)A + B{t) is the Unear operator studied in Section H] and 

F{t,w,Aw) = {^'{v4t) + wy^ -^'{v^{t))-^)D{^'{v4t)+w)) Aw 

+ ^'iv4t))-^D {^'(v^t) +w)- ^'Mt))) Aw , 
G{t,w,Vw) = ^'{v^t) + wy^D'i'"{v4t) + w)[Vw,Vw] , 

H{t,w) = {0,g2it,w))'^ , where g2{t,w) = g2{v^{t) + w)[wi„Wh] . 

Clearly, F{t,w, Aw) = 0{wAw) and G{t,w,Vw) = 0{\Vw\'^). More precisely, if ||?y||L°o is suffi- 
ciently small, we have 

\\F{t,w,Aw)\\Li < C\\w\\l2\\Aw\\l2 , \\G{t,w,Vw)\\Li < C\\Vw\\l, , 

\\F{t,w,Aw)\\L2 < G\\w\\l^\\Aw\\l2 , \\G{t,wyw)\\L2 < C\\Vw\\l,. 



Since 
we find 



|Vi(;|||2 < C||w||j^2 II A7i;||j^2 and ||Vit;||^4 < C||t(;||Lcx) || At(;||^2 



(5.2) 



\\F{t,w,Aw)\\Li + \\G{t,w,Vw)\\Li < G\\w\\l2\\Aw\\l2 , 
\\F{t,w,Aw)\\L2 + \\G{t,w,Vw)\\L2 < C7||w;||i^||A'u;||i2 . 
Under the same assumptions, we also have 

\\H{t,w)\\Li < C||tt;h||^2 and \\H{t,w)\\L2 < C||?x;h||L2 ll^^hlU"" • (5.3) 

Setting K(t, w, Vw, Aw) = F{t, w, Aw) + G{t, w, Vw), we can write the integral equation associated 
with (|3.6p in the form, 

w{t) = M{t,0)wo + [ M{t,s)K{s,w,Vw,Aw)ds+ [ M{t, s)H{s,w) ds . (5.4) 
Jo Jo 

We now describe the function space in which we shall look for solutions of (j5.4p . Assume that 
n < 3 and choose r G (4, +oo), so that 

1 n 1 , , 

- < - < 1--. 5.5 
r 4 r 

We define the Banach space 

Y = {ti;eC°([0,+oo),Li(M")nL^(M"))nL^X(0,+oo),ij2(Mn))| ||^„||^ < oo} , 

where 

||w;||y = sup||w;(t)||ii +sup(l + t)"/2||u;(t)||^^ + / (1 + t)'^||Au;(t)||^2 dt . (5.6) 
t>o t>o V^o / 
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Lemma 5.1 If wq G L^(]R") n H'^{W^), the linear solution Wo{t) = M{t,0)wo belongs to Y and 
\\Wo\\y < Ci(||wo||li + ll'^^oll/fs) for some Ci > 0. 

Proof. Since n < 3, we have H'^{R"-) ^ L°°{W), hence wq G ^^(IR") n L~(M"). Thus, using the 
hnear estimates estabhshed in Proposition 14.41 obtain 

snp\\Wo{t)\\Li+snp{l + tr/^Wo{t)\\L^ < C(\\wo\\li + \\wo\\l^) . 

t>0 t>0 



On the other hand, using Proposition 14.41 and Remark 14.51 we find 

'tl+n/4 



II II 

AT^oWIIl^ < C\\w^\\n^ fort<l, ||AWo(i)|U2 < '["^^ f for t > 1 



As rn/4 > 1 by (jS.Sp . we conclude that 

l/r 



(i+tniAw^o(t)irL2dtj < c{\\w^u.^\\w^\\u^). 



The next step consists in estimating the integral terms in (15. 4p . namely 

Z{t) = [ M{t,s)K{s,w,Vw, Aw) ds , and J{t) = [ M{t, s)H{s,w) ds . 
Jo Jo 

Proposition 5.2 There exist C2 > and 62 > such that, for all w £Y with \\w\\y < S2, we have 

llXlly + IUIIy < Calklly. 

Proof. Throughout the proof C denotes a constant that changes between estimates, but does not 
depend on w. Smallness of in y implies that estimates (j5.2p and (j5.3p hold for the nonlinear ities. 

Estimate on ||J(t)||j;^i 

Using ()4.19p with p = q = 1 and the first estimate in ()5.2p , we find 

\\I{t)\\Li<C ['\\w{s)\\L2\\Aw{s)U2ds 



Jo 

< C\H\y f 7-^^(1 + s)\\Aw{s)\\l^ ds (5.7) 

In the second inequality we used the bound II 1 1 2,2 < ||'w(s)||]/i^||w^(s)||/^(^ < ||ii;||y(l + s)~"/^, and 
in the last line Holder's inequality. Taking the supremum over t > and using (j5.6p . we conclude 
that 

sup||T(t)||ii < C\\w\\l. 
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Estimate on ||T(t)||Loo 

For t < 1, we use (|4.19|) with (p, q) = (2, oo) and the second estimate in (|5.2p . We obtain 



ITOIU^ <CyJ ^^_^^^^^^J w{s)\\L^\\Aw{s)h2 ds<Ct'-'i-^\\wrY, (5.8) 
where the last estimate is again a consequence of Holder's inequality. Note that 1 — — ■^>Oby 



(j53]) . For t > 1 we split 

l-t/2 ft 

l{t)= I M{t,s)K{s,w,Vw,Aw)ds+ / M{t, s)K{s,w,Vw, Aw) ds 

Jo Jt/2 

= : Ii{t)+l2{t). 

Using (I4.19P with {p, q) = (1, oo) and proceeding as in ()5.7p . we find 

(1 + tTl'\\Ut)\\L^ < C{1 + tr'^ __i_||^,(,)||^, \\^W{S)U. ds 

l-t/2 

<C \\w{s)\\L2\\Aw{s)\\L2ds <C\\w\\\ . (5.9) 







On the other hand, using ()4.19p with (p, q) = (2, oo) we obtain 

{l + tT''\\l2{t)\\L-<C{l+tr" f —\j-^\\w{s)\\L^\\Aw{s)U.ds 

Jt/2 [t - S)"/* 

- ^"""^ (TTT) i/l (T^^^' + s)\\Awis)U. ds (5.10) 



<C\\wfy- --r, 

(l+t)4+f 

where in the last line we used Holder's inequality as in (j5.8p . This shows that 

sup(l + t)"/2||j(t)||^^ < C7||u;||^. 

t>o 

Estimate on ||AT(t) 11^,2 

The estimates for the second derivative require maximal regularity, Proposition 14. 7i We need to 
estimate /o°°(l + i)'^||^2r(t)||^2 dt. We therefore split the integral and first estimate 

ril + tY\\AIit)\\l,dt < C r \\witWL^\\Awit)\\l,dt<C\\wfY\ 
Jo Jo 

where we used Proposition 14.71 and estimate (15. 2p in the first inequality, and the definition of ||y 
in the second inequality. 

We next derive estimates for AI{t) at t >T. Here it is more convenient to decompose 
I(t) = [ M{t,s)K{s,w,Vw,Aw)ds+ [ M{t, s)K{s,w,Vw, Aw) ds 

Jo Jt-T 

=: I3{t)+l4{t). 
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Using (I4.19P with (p, q) = (1, 2) and Remark 14.51 we control the first term as 

{l + t)\\AUt)h2 < C{l + t) f ^ ^ \\w{s)\\L2\\/\w{s)\\L2ds 

Jo [t-sj'-^i 



t-T ^ 1 



< C\\w\\y{1 + t) / ^(1 + s)\\Aw{s)h2 ds 

Jo [t - S)^^4 (1 + s)^^4 



„ „2 1 



T 



a+t) 

where in the last estimate we used Holder's inequality together with the fact that 

f /"V 1 1 1 \« , c , , , 

U<^i + t-.(t-.)"/4(i + ,)i.fj ^(TTT)^' ^'-''^ 

Integrating over time and recalling that nr/4 > 1, we obtain the desired bound for X3: 

{i + ty\\AUt)\\i,dt < C\\wfy\ 

The other term 24 is estimated directly using (j5.2p and Corollary I4.8t 

/ ii+trwAumi.dt < c / {i+tY\\w{t)\\i^\\Awmhdt < cwwfy'. 

JT Jo 

Estimates on \\J{t)\\Li and \\J{t)\\Loo 

In this term, the nonlinearity does not contain derivatives which would yield decay, but we can 
exploit the the stronger decay of the linear operator J^{t,s) when acting on H{s,w). Indeed, since 
H{s,w) = {0, g2{s,w))^ , we have 

\ Mhh{t,s)g2{s,w) J 
and it follows from Proposition 14.41 that 

\\M{t,s)H{s,w)h,iU'^) < C—^ u;)||iP(Rn) , (5.12) 

for 1 < p < q < 00. Using (|5.12|) with p = q = 1 and estimate (|5.3p . we thus find 

IIJWIIli < C f \\w{s)\\l,ds < C\\wfy f \ ds < C\\wfy. (5.13) 

jQl + t-S Jq 1 + t - s (1 + s)"/^ 

In a similar way, using (j5.12p with {p,q) = (2,oo), we arrive at 

(i + (r/'||j(i)llL» < C(i + 1)"/" £ 1 ^ i _ , (t _',),./4 \M')h~\\'^(')\\L' tia 

< + tr'-' jl (5.14) 



< C7||u;ii2 



Y 
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Estimate on ||AJ'(t)||i2 
We first observe that 



AJ(t) = / M.h{t,s)Ag2{s,w)ds. 
Jo 



Since g2{t, w) = g2{v^{t) + ■w)[wii, tt^h], it is clear that \Ag2{t, w)\ < C(|i(;| | Aw| + | Vwp), so that Ag2 
satisfies the same estimates (|5.2p as F and G. Using the L^-L^ estimate for Ai.h{t, s), we thus find 

< C\\w\\l- 

- II 11^^(1 + i)n/4 ' 

where in the last line we used Holder's inequality and estimate (jS.lip . We deduce that 

/•oo 

/ {l+tr\\AJ{t)\\l,dt < C\\wfy\ 

Jo 

and the proof of Proposition 15.21 is complete. ■ 

Proof of Theorem [1] (n < 3). As was observed in Section [3l we can work with the transformed 
equation (j3.6|) instead of the original perturbation equation (j2.ip . Also, we can assume without loss 
of generality that the initial perturbation wq satisfies H^yollii + ||^^o||/f2 ^ <^o for some small 6o > 0. 
Under these assumptions, we can solve equation (15. 4p by a standard fixed point argument in the 
Banach space Y defined by (15. 6p . Indeed, let Af denote the right-hand side of ()5.4p . namely 

{Afw){t) = Ji4{t,0)wo+ [ J^{t,s)K{s,w,Vw,Aw)ds+ [ J^{t, s)H{s,w) ds . 

Jo Jo 

w gY satisfies \\w\\y < ^2, where (52 > is as in Proposition 15.21 we know that Mw E Y and that 

\\Mw\\y < Ci6o + C2\\w\\lr . (5.15) 

Similar calculations show that 

\\Afw-Afw\\Y < C2{\\w\\y + \\w\\y)\\w - w\\y , (5.16) 

whenever w,w G Y with ||w||y < ^2, H^iHy < ^2- Let B G Y denotes the ball of radius R = 
min(2Ci(5o, (^2) centered at the origin. If 5o > is small enough so that 2C2R < 1, it follows easily 
from (j5.15p . (j5.16p that N{B) C B and that is a strict contraction in B. Let w gY he, the unique 
fixed point of N in B. Then w is a global solution of ()3.6p . and if we return to the original variables 
by setting u{t,x) = ^'(f*(t) + 'w(t,x)), we obtain a global solution of (11. ip which satisfies the decay 
estimate ()1.6p (with to = 0), because 

\\w\\y 

\uit, x) - u.{t)\ = l^iv^t) + wit, x)) - ^iv,it))\ < C ^^^ l^^i2 ' ^ > • 
This concludes the proof. □ 
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Remark 5.3 The limitation n < 3 in the above proof is due to the fact that we use maximal 
regularity (MR) in the Hilbert space L^(R") only, see e.g. ()5.8p . This choice was made for simplicity, 
but for equation (j4.2p it is known that MR holds in all spaces with 1 < p < oo, see \19^ . It 
is not difficult to verify that the argument above can be adapted to any space dimension n if we use 
MR in L^iW^) with p sufficiently large, depending on n. 



6 Asymptotic behavior 

We know from Theorem [1] that small, localized perturbations of the periodic solution n*(t) converge 
to zero like as t — > +oo. This decay rate is optimal in general, and it is even possible to compute 
the leading term in the asymptotic expansion of the perturbation as t — > +oo. In this section, we 
assume (for simplicity) that 1 < re < 3 and we consider the solution w{t, x) of (13. 6[) with small 
initial data wq S //-^(M") niJ^(M"'). If we decompose this solution as w{t,x) = {wQ{t,x),w\i{t,x)Y' , 
we first observe that the hyperbolic part w\i{t,x) € M^^"*^ has a faster decay as t ^ oo. 

Proposition 6.1 If the initial data wq E L^(M") n i?^(M") are sufficiently small, the hyperbolic 
component of the solution w of (j3.6p satisfies 



sup(l + ^)||^i;h(^)||Ll +sup(l + t)i+?||?i;h(OllL- < <^(lko||Li + Ikolk^) • 

t>0 t>0 
Proof. Projecting the integral equation (j5.4p onto the hyperbolic component, we find 

w\,{t) = M\,.{t,{))wQ+ I M\,.{t,s)K{s,w,Vw,/^w)ds + I Mh-{t, s)H{s,w) ds 

Jo Jo 

=: MUt,0)wo +Ii,{t) + Jh{t) , 

where A^h-(*;'S) = {Mhoit, s), Mhh{t, s)). We know from Proposition 14.41 that 

1 1 

\\Mh-{t,s)w\\Li(Rn^ < C — — n,i 1, \\w\\lp{R") , (6-1) 

i + t - S _ g^aVp q) 

for 1 < p < q < oo. In particular, we have 

sup{l + t)\\Mh-it,0)wo\\Li + sup{l + t)^+^\\Mh.{t,0)wo\\L^ < Cdlu-ollii + H^i'oIIl-)- 

t>0 t>0 

Moreover, proceeding as in the proof of Proposition 15.21 we obtain 

||Xh(t)||Li < C f—^ "^''^ {I + s)\\Aw{s)U. ds <^ \\w\\l, 

Jo 1 + t- S (1 + s)^+4 1 + t 

and the same result holds for (1 + t)"/^||Th(i)||L°o. Finally, to bound the term Jh, we observe that 
Aih-it^ ■s)-ff (s, w) = M.hh{t, s)g2{s, w) and we use the strong decay in time given by Proposition 14. 4[ 
We thus find 

||Jh(t)||Li < C / — -.\\wi,{s)\\Li\\wi,{s)\\L^ds (6.2) 



(l + t-s)2 

/ ^ /■* 1 ll^lly ^ / ^11 112 

- Jo (1 + ^-5)2(1 + 5)^^/2 -(l+i)n/2ll 
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and the same result holds for (1 + t)"/2||Jh(t)||L°°- This gives the desired result if n > 2. If n = 1, 
we only have 

sup{l + t)^/^\\wh{t)\\Li +sup{l + t)\\wh{t)\\L^ < CWwoh^nH^ , 

t>0 t>0 

but if we now return to (j6.2p we obtain the stronger estimate 

which also holds for (1 + t)^^^||,7h(i)||L°°- This concludes the proof. ■ 

We next consider the central component WQ^tjx) G M, and prove that it behaves asymptotically 
like a solution of a linear equation with suitably modified initial data. 

Proposition 6.2 If the initial data Wq G L^(W^) n H'^(W^) are sufficiently small, the central com- 
ponent of the solution w of (13.61) satisfies 



(J 

\\wo{t) - MoQ{t,Qi)woo\\L^ + (1 + t)"/^||u7o(t) - MoQ{t, 0)ti;oo||L- < (Y^^ (IkollLi + WmWm) , 
where 7 = | + i < 1 and w 00 e L1(M") n L°°(M") is defined by 



Woo = {^0)0+ / Ko{s,w{s),Vw{s),Aw{s))ds. 
Jo 

Proof. Projecting (|5.4|) onto the central component, we find 

wo{t) = Mo.{t,0)wQ+ f Mo.{t,s)K{s,w,Vw,Aw)ds+ I Mo.{t, s)H{s,w) ds 

Jo Jo 

=: Mo.{t,0)wo+Io{t) + Jo{t), 

where M.o.{t, s) = {M.oo{t, s),M.ohit, s)). Our goal is to extract from wo{t) the leading contributions 
as t — > +00. The last term ,Jo{t) is clearly negligible in this limit. Indeed, using Proposition 16.11 and 
proceeding as in (|5.13p . (|5.14p . we obtain 

\\Jo{t)U, + {l + tr/^Mt)\\L^ < ^^Wwoh^nH^. (6.3) 

The same estimate holds for the linear term A^oh(i) 0)(wo)h) because A^ohC^iO) decays as fast as 
(1 + t)~^A4oo{t,0) by Proposition 14.41 Using the same remark and proceeding as in (j5.7p . (j5.9p . 
(jS.lOp . we see that the integral term Xoh(i) := Jq M.oh(t, s)Kh{s,w,Vw, Aw) ds also satisfies ()6.3p . 
So the only remaining terms are M.oo{t,0){wo)Q and 

^oo{t) = / Moo{t,s)Ko{s,w,Vw,Aw)ds 
Jo 

ft ft/2 
= / Moo{t,s)Ko{s,w,Vw,Aw)ds + / {Moo{t, s) - Moo{t,0))Ko{s,w,Vw, Aw) ds 

Jt/2 Jo 



+ A^oo(i,0) / KQ{s,w,Vw,Aw)ds - A^oo(i,0) / Ko{s,w,Vw,Aw)ds 

Jo Jt/2 

--: Ioi{t) +Io2{t) +loz{t) +loi{t) ■ 
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Proceeding as in ()5.7p . (IS.lOp . it is straightforward to verify that 

iiJoi(t)iiii + (i+tr/2||jo,(t)||^^ < _^_||^||^, t>o, 

where 7 = ^ + ^ < 1, and the same estimate clearly holds for Xo4(i) too. Finally, using Remark 14.61 



to bound the difference A4oQ{t, s) — M.oo{t, 0), we obtain 

||Xo2(t)||Li + (l + t)"/'||Xo2(t)||L- < C 1^^^ j\\Ko{s,w,Vw,Aw)\\Lids < j^^^Wwfy. 



This concludes the proof of Proposition 16.21 because Aioo{t, 0)(tt;o)o +^03{t) = -^00 (^) 0)woo- ■ 

It is now rather easy to prove Theorem [2] (in the case where n < 3) . Combining Propositions 16.11 
and 16.21 we find 

\\wit) - eiWmL^ + a + tr^^wit) - eiWit)\\L^ < .^-^(\\ujohi + \\nJo\\H^), (6.4) 

where W{t) = M-Qo{t,0)woc, and ei = (1,0)-^ is the first vector of the canonical basis in M^. 
Furthermore, we claim that 

||r/2p7(i,xt^/2) -dGlliinioo > 0, (6.5) 

where G is defined in (jl.7p and 

a= WoQ{x)dx = { WQ{x)dx+ / / K{t,w,'Vw,Aw)dxdt\ . (6.6) 

To prove the L°° claim in (16.51) . we use Fourier transforms and simply note that the quantity 

\\W{t,kt-'/^)-dG{k)Ui = ||Moo(t,0;A;t-i/2)tI;^(A;t-i/2) -e-^o'^'u;oo(0)||Li 

converges to zero as t ^ 00 by Lebesgue's dominated convergence theorem, in view of Proposition ITT] 
and Remark 14.21 The claim can be established in a similar way, using the same ideas as in the 
proof of Proposition 14.41 (we omit the details). 

We now return to the original variables. Since M^=(t) = ^{ujtei), the solution of (jl.ip given by 
u(t,x) = 'if{v(t,x)) = ^(z;^,(t) + w{t,x)) can be decomposed as in (jl.Sp . with a{t,x) = uj~^W{t, x) 
and 

p{t, x) = ^'(v=,(t) + w{t, x)) - ^'(v*(t)) - 'if'{v4t))w{t, x) + ^''(w*(t))(u;(t, x) - eiW{t, x)) . 

Estimates (|6.4p . (j6.5p immediately give (jl.9p . (jl.lOp with = uj~^a. Finally, the formula (jl.lip for 
a=K follows from the expression (16. 6p of d and the fact that ?7*(0) = (^'(0)~^)^ei, u^O) = a;^'(0)ei. 
This concludes the proof of Theorem [2l □ 



7 Examples and perspectives 

In this final section, we first give a simple example of a 2-species reaction-diffusion system with 
a periodic orbit n*(t) which is asymptotically stable for the ODE dynamics but does not satisfy 
Hypothesis II. 2[ We then discuss possible generalizations of the results of this paper. 
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7.1 Destabilization by diffusion: a simple example 

One may feel inclined to believe that ODE-stable periodic orbits tend to be stable for the PDE 
dynamics, that is, that our Hypothesis 11.21 is satisfied in most cases where the periodic orbit is 
stable for the ODE. Our example below shows that this is not the case, even for a simple reaction- 
diffusion system with only two species. 

We consider the following 2-species reaction-diffusion system 

ut = Duxx + Ju + (e^ — \u\^)Ru , (7-1) 

where u = {ui,U2)'^ £ ffi^ and |ttp = uf + u\. Here e > is a parameter, D is a 2 x 2 real matrix 
with positive eigenvalues, and 

/^O -A o lcos{e) -sin(^)A ^ , ,^ 

'=[l oj' ^=(^smW cos(Vj' 

The system (|7.ip has a 27r-periodic solution = eu{t), where u{t) = (cos(t), sin(t))^. Linearizing 
(j7.1|) at u^{t) we obtain 

vt = Dvxx + iJ-2e^Ru{t)u{tf)v, 

or equivalently 

Vt = -k^Dv + ( J - 2e^R u{t) u{tY)v , A; G R . (7.2) 

Of course v{t) = u'{t) is a solution of (j7.2l) for k = Q. 

Let F = Fo(27r, 0) be the Floquet matrix associated to (|7.2p for A; = 0. Then 

Det(F) = exp(^Tr j {J -2e^Ru{t)u{tf)dtj = exp(-27re^Tr(i?)) . 

The Floquet exponents (for = 0) are therefore Ai = and A2 = — e^Tr(i?). As Tr(i?) = 2cos(0) > 
0, it follows that n^,(t) is a stable periodic orbit for the ODE dynamics associated to ()7.ip . 

To compute the Floquet exponents for small k, we consider the adjoint ODE 

Ut = {J + 2e^u{t)u{tf R^)U . (7.3) 

As is easily verified, the unique nontrivial bounded solution of (j7.3p is f7*(t) = Ru'{t). Using 
formula p.4p . we conclude that 

If the diffusion matrix D is symmetric, then Tr(L') > and Tr( JD) = 0, hence necessarily do > 0, 
which means that the periodic solution u^{t) is spectrally stable for long wave-length perturbations. 
But if Z) is a nonsymmetric matrix, then Tr{JD) 7^ and therefore we can choose 6 G (— 7r/2,7r/2) 
in such a way that do < 0. This gives an example of a periodic orbit exhibiting a sideband instability. 
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On the other hand, as e ^ 0, the periodic orbit reduces to the fixed point u = 0, for which 
it is easy to perform a stabihty analysis. For a fixed wavenumber fc G M, we have to compute the 
eigenvalues Ai(A;), X2{k) of the linearized operator J — Dk^. By direct calculation, we find 

Tr(J-DA;2) = -Tr(D)A;2 < 0, and Det{J - Dk"^) = 1 + TT{JD)k'^ + Det{D)k^ . 

If we choose D such that Tr(JL>) + 2(Det(L>))V2 < 

we see that there exists a nonempty open 
interval / C (0, +oo) such that Tr( J - DP) < and Det( J - DP) < if k^ £ I. Thus, one 
of the eigenvalues Xi{k) is strictly positive, which means that the equilibrium ti = is unstable 
with respect to perturbations with wavenumbers k such that k"^ S /. By continuity, this Turing 
instability persists for the periodic orbit n*(t) = eu{t) if e > is sufficiently small: for k'^ G /, one 
of the Floquet exponents has positive real part. 

Summarizing, the periodic solution = eu{t) of the reaction-diffusion system (jT.ip exhibits: 

(i) a sideband instability, if Tr(L') — tan(0)Tr( JD) < 0; 

(ii) a Turing instability, if Tr(Ji:>) + 2(Det(i:'))^/^ < and e < 1. 

Remarks 7.1 

1. The instability criteria above are never satisfied if the matrix D is symmetric. But if D has 
eigenvalues di > d2 > 0, we can choose an invertible matrix S so that S^^DS = V = diag((ii, ^2). 
Then setting u = Sw we obtain the equivalent system 

wt = Vw^^ + S-^JSw + S-^{1 - \Swf)RSw , (7.4) 

where now the diffusion matrix has the usual, diagonal form. The Floquet exponents characterizing 
the stability properties of the periodic orbit are of course unaffected by this linear transformation. 
Thus we can find 2- species systems of the form (17. 4|) . with diagonal diffusion matrix, which exhibit 
either a sideband or a Turing instability. 

2. Our example is clearly reminiscent of the complex Ginzburg-Landau equation ( CGLE), 

ut = (1 + ia)An + ti — (1 — «c)|iipu , (7.5) 

where a,c are real parameters and u : x C. The system ()7.5p possesses a homoge- 

neous time-periodic solution of the form u{t,x) = e*'^*, which exhibits a sideband instability if ac > 1 
(T3enjamin-Feir criterion^ and a Turing instability in other parameter regions, see e.g. [1]. The com- 
plex Ginzburg-Landau equation arises as a modulation equation near Hopf bifurcations in reaction- 
diffusion systems (see, for example, l2T\ [T7|/j. One therefore expects stability and instability proper- 
ties of small- amplitude periodic solutions near Hopf bifurcation to be governed by those of the CGLE; 
see, for example, 122] for a result in this direction. Since the diffusion matrix of CGLE possesses 
complex eigenvalues, it cannot be cast as a real reaction- diffusion system with diagonal diffusion 
matrix. Our example and the remark above show that one can almost explicitly recover the proper- 
ties of CGLE with diagonal diffusion matrices. Upon substituting the diffusion matrix of CGLE in 
our example, one would recover precisely the Benjamin-Feir criterion from the instability criterion 
Tr{D) — tan{6)Tr{JD) < 0. In a different direction, one could also extend our results to diffusion 
matrices with D + > without any additional difficulties, which would then include CGLE as 
a particular example. In the context of reaction- diffusion modeling, cross- diffusion phenomena that 
are associated with off-diagonal elements of D are however quite uncommon. 



26 



7.2 Discussion and perspectives 

We believe that the method presented here can be adapted to other situations. We mention the 
stabiUty of wave trains, u{kx — ujt), with u{(,) = u{(^ + 2tt) and uj,k > 0, and Turing patterns 
u{kx) = u{—kx), with u{S^) = + 27r) and A; > 0. In both cases, one finds continuous spectrum 
with diffusive decay properties for the Unearization. In both cases, the absence of relevant, self- 
coupling terms in the neutral mode has been shown previously; see [2b\ [S] . 

Interesting questions arise when one attempts to extend the class of allowed perturbations. One 
may for instance consider perturbations v such that Vv € L^. In one space-dimension, this would 
correspond to perturbations with different phase shifts at x = ±00. One would still expect the 
diffusive linear part to be dominant so that one would find error function asymptotics for the phase 
correction. 

In fact, one would expect some type of stability for much more general perturbations. For 
instance, in one space-dimension, the homogeneous oscillation u^(t) is embedded in a family of 
wave train solutions u{kx — ujt;k), with ~ and u ~ u^; [25 1 Section 3.3]. Under our stability 
assumptions, using Lyapunov-Schmidt reduction, one finds m = io{k) = loq + uj2k'^ + 0{k'^); see for 
instance j20|, Lemma 2.1]. Of course, u{kx — ut; k) is not close to any fixed homogeneous oscillation, 
but it is close to an appropriate phase shift of the oscillation in any finite region of space. All solutions 
in the basin of attraction of such wave trains then stay close to our homogeneous oscillation, orbitally, 
and pointwise in space. A natural question then asks for the asymptotics of initial conditions of the 
type u{k{x)x — uj{k{x))t;k{x)), where k{x) — > k± for x ±00. In the case of the real Ginzburg- 
Landau equation, which exhibits spatial oscillations u{kx), this question was addressed in [6l [TT]. 
showing that asymptotics are governed by a nonlinear diffusion equation 9t = W{9x)x, with locally 
uniform convergence to a fixed, intermediate wavenumber. Near temporal oscillations, one expects 
dynamics to be governed by a viscous conservation law 6t = d{6x)x + j{Gx), so that solutions with 
asymptotically constant wavenumber would be expected to converge to viscous shocks [HI Section 
8], or to rarefaction waves. 

More generally, one could ask about the orbital stability of a family of oscillations: starting with 
an initial condition u{k{x)x — Lo{k{x))t; k(x)), with k{x) bounded and k' small, will the solutions 
remain close to a solution of that form for all times? Again, this question possesses a simple 
answer for the approximation of the dynamics by a viscous conservation law due to the maximum 
principle, which gives immediate supremum bounds on 6x in terms of the initial condition. For 
the full reaction-diffusion dynamics however, which are only approximately described by a viscous 
conservation law, this question remains wide open. 
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